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Recently F. Huang [Commun. Theor. Phys. 42 (2004) 903] and X. Tang and P.K. Shukla 
[Commun. Theor. Phys. 49 (2008) 229] investigated symmetry properties of the barotropic 
potential vorticity equation without forcing and dissipation on the beta-plane. This equation is 
governed by two dimensionless parameters, F and /3, representing the ratio of the characteristic 
length scale to the Rossby radius of deformation and the variation of earth' angular rotation, 
respectively. In the present paper it is shown that in the case F ^ there exists a well-defined point 
transformation to set (3 — 0. The classification of one- and two-dimensional Lie subalgebras of the 
Lie symmetry algebra of the potential vorticity equation is given for the parameter combination 
F ^ and f3 = 0. Based upon this classification, distinct classes of group- invariant solutions is 
obtained and extended to the case /? 7^ 0. 
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1 Introduction 



2 The Lie symmetries 



There is a long history in dynamic meteorology to 
use simplified models of the complete set of hydro- 
thermodynamical equations to gain insides in the 
different processes characterising the various struc- 
tures and pattern occurring in the atmosphere. One 
of the most classical models in atmospheric sci- 
ence is the barotropic (potential) vorticity equa- 
tion. It has been successfully used both for theo- 
retical considerations [181 [7] and practical numer- 
ical weather predictions |S] since it is capable of 
describing some prominent features of mid-latitude 
weather phenomena such as the well-known Rossby 
waves and blocking regimes. In nondimensional 
form it reads [T7] 
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Here "0 is the stream function, C, the vorticity, 
J{a,b) ~ da/dxdb/dy — da/dydb/dx is the Jaco- 
bian and F and j3 represent the ratio of the charac- 
teristic length scale to the Rossby radius of defor- 
mation and a parameter describing the variation of 
earth' angular rotation, respectively. 



The symmetries of ([T]) have first been investigated 
in [S] and were applied to construct new solutions 
from known ones. Later, they have been studied 
by [To] and [H] , but without reference to the group 
classification problem. This problem arises since 
different values of F and (3 lead to different sym- 
metry properties of ([T]), which in turn characterise 
different physical properties of model H]). However, 
as shown in this paper, there are only three essential 
combinations of the values of these two parameters, 
given by F = 0, /3 = 0; F = 0, /3 7^ and F 7^ 0. 
The first combination leads to the usual vorticity 
form of Euler's equation which is the issue of e.g. 
[2]. This combination is of particular interest, since 
it gives rise to a new symmetry [3] (a so-called po- 
tential symmetry) which is not present in the veloc- 
ity form of Euler's equation. The second combina- 
tion of parameters was discussed in [H [TT] . It is the 
usual barotropic vorticity equation. The parame- 
ter (3 can be set to 1 by scaling and/or changing 
signs of variables. Note that for both combinations 
the associated Lie invariance algebras are infinite 
dimensional. 

Now we show that if F 7^ 0, we can always set 
(3 — Q. (Then the nonvanishing parameter F can be 
scaled to ±1.) For this purpose, we recompute the 



symmetries of ((!]) for the case F ^ Q and (3 arbi- 
trary. This is done upon using the computer algebra 
packages MuLie [S] and DESOLV 0. For this com- 
bination, equation ([l]) admits the six-dimensional 
Lie symmetry algebra generated by the opera- 
tors 



V = tdt- i-td. 
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This algebra is not singular in (3 and consequently 
it also includes the case /3 = 0. Moreover, comput- 
ing the symmetries for the case /? = explicitly, we 
obtain the same algebra (U) with /? = 0. Hence 

V = tdt~ i^d^, 
Vr = -2/<9x + xdy, 
Vt = dt, 
Vx = dx, 

Wy ^ dy, 

— dif, 

are the generators of the Lie symmetry algebra Qq. 
The physical importance of these generators is the 
following: V generates simultaneous scaling in t and 
"0, Vj, is the rotation operator in the (x, y)-plane, 
Vt, Vj,, Vy and v^ are the infinitesimal generators 
of translations in t, x, y and ip, respectively. The 
nonzero commutation relations between basis ele- 
ments ([3]) are exhausted by 

[vt,23] Vt, [v,0,2?] = -v^, 

[Vx,Vr] Vy, [Vy,Vr] = -Vx- 

Therefore, the algebra Oq has a simple structure. It 
is a solvable Lie algebra and can be represented as 
the direct sum © e(2), where e(2) = (v-c, v.y, v^} 
is the Euclidean algebra in the (x, ?/)-plane and 
03 4 = (vt,v^,2?) is a three-dimensional almost 
Abelian Lie algebra from Mubarakzyanov's classifi- 
cation of low dimensional Lie algebras [13] . 

It is straightforward to show that the Lie algebra 
maps to the Lie algebra ^ under the transfor- 
mation given by 



f3 

t = t, x^x + —t, y = y, 
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This transformation also maps ([T]) to the same equa- 
tion with /3 0. That is, (g]) is an equivalence trans- 
formation for the class of equations of the form ([1]) 



with F ^ 0. Hence the /3-term can be neglected un- 
der symmetry analysis. Every solution of the equa- 
tion with /3 = can be extended to a solution with 
/? 7^ by means of the transformation ([4]). 

We also note that the maximal (infinite dimen- 
sional) Lie symmetry algebras in the cases F = 0, 
/? = and F = 0, /3 7^ are neither isomorphic 
to each other nor isomorphic to the (finite dimen- 
sional) algebra Qq. Consequently, it is not possible 
to find point transformations that relate the corre- 
sponding PDEs to each other. All the transforma- 
tions used for the reductions of the parameters F 
and (3 belong to the equivalence group of class ([T|). 
The group classification list for class ^ is therefore 
exhausted by the three inequivalent cases F = 0, 
/3 = 0; F = 0, /3 = 1 and F = ±1, /3 = 0. 

3 Classification of subalgebras 

Classification of subgroups of Lie symmetry groups 
of differential equations is an essential part in the 
study of these equations. This is since classifica- 
tion allows for an efficient computation of group- 
invariant solutions, without the possibility of an oc- 
currence of equivalent solutions. Classifying sub- 
groups may further lead to the construction of 
simple ansatze for the corresponding equivalence 
classes of reduced differential equations. Thereby, 
the classification also provides an important step for 
further investigations of properties of these reduced 
equations. 

The classification of subgroups of symmetry 
groups is usually done by the classification of the 
associated Lie subalgebras with respect to the ad- 
joint representation [Ml [15]. The potential vorticity 
equation ([T]) is a (2-1-1) model and thus Lie reduc- 
tions up to an ordinary differential equations re- 
quire the classification of one- and two-dimensional 
subalgebras. An exhaustive classification of subal- 
gebras exists only for Lie algebras up to dimension 
four [16]. Thus we need to classify subalgebras of 
Qo . This problem is not difficult because the algebra 
ao has a simple solvable structure. 

The adjoint representation of a Lie group on it's 
Lie algebra is given as the Lie series 



w(e) = Ad(e"^)wo 



oo „ 

y- 

Tl = 



{v",wo}, (5) 



with {•, •} being defined recursively: 
{v°. Wo} := Wo, 

{v",wo} :=(-l)"[v,{v"'i,wo}]. 

Alternatively, the adjoint representation can also be 
calculated by integrating the initial value problem 

dw(e) 



d£ 



[w(e),v], w(0) = wo. 
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For basis elements ([3]) of the algebra oo we obtain 
the following nontrivial adjoint actions 



Ad(e"^')P 




Ad(e'^^'")P 


= V + ev^ , 


Ad(e^^)vt 




Ad(e^^)v^ 




Ad(e-'^-)v^ 


= V,. - evy, 


Ad(e^^«)v,. 


= V,. + eVa;, 


Ad(e-'^'-)v, 


= Va; cos £ + Vj, sin £ 


Ad(e-'-)v, 


= — Va; sin £ + Vj, COS £ 



3.1 One-dimensional subalgebras 

The classification of one-dimensional subalgebras of 
the whole symmetry algebra ^ is done by an induc- 
tive approach We start with the most general 
infinitesimal generator, 

V = aiV + a2V,. + 03 vt + 04 v^, + a^Vy + aev^, 

and simplify it as much as possible by means of ad- 
joint actions. Depending on the respective values of 
the coefficients a;, i — 1, . . . , 6, we find the follow- 
ing list of inequivalent one-dimensional subalgebras 
of @: 

1. ai ^0,a2^ 0: 

(V + aVr). 

2. ai 7^ 0,02 = 0,04 ^ 0: 

{V + av,). 

3. ai = 0, 02 7^ 0, 03 ^ 0: 

(v,. ± vt + av^,). 

4. ai = 03 = 0, a2 7^ 0: 

(6) 

(Vr + cv^). 

5. fli = 02 = 0, 03 7^ 0, 04 7^ : 

(vt + avj, + cv^,). 

6. ai = 02 = 03 = 0, 04 7^ : 

(v^, + cv^,). 

7. ai = a2 = a3 = 04 = 05 = : 

where c G { — 1, 0, 1} and a, G M. In case 5 we can 
additionally set a G { — 1,0, 1} if c = 0. 

3.2 Two-dimensional subalgabras 

The classification procedure of two-dimensional 
subalgebras follows in essential the same way as 



the one-dimensional case. The two most general 
linearly independent infinitesimal generators 

ry r\ r\ ry r\ r\ 

V = a^V + 02 + 03 Vt + a^Vx + a^Vy + a^-v^ 

are simultaneously subjected to the adjoint actions 
and nonsingular linear combining under some as- 
sumptions on the coefficients , i = 1, . . . , 6, j = 
1,2. Moreover, the required closure property of the 
subalgebra (v^,v^) with respect to the Lie bracket 
(i.e. [v^,v^] G (v^,v^)) eventually places further re- 
strictions on the coefficients. By applying this tech- 
nique, we find a set of inequivalent two-dimensional 
subalgebras of ([3]). For reason of brevity, we only 
list the subalgebras without the corresponding con- 
ditions on the respective coefficients a\ : 

(2?,v,.), 
{V + avr, Vf), 
{V + avj,, Vi), 
{V + av.j;,Vy), 
{V + avr, v^,), 

(P + aVa^V^), 
(Vr + CVt, V^}, 

(vt + av^ + cv^,,^fy + 6v^), 

(vt -I- av.^, v^), 

{■v.j, + cv,i,,\-y + 6v^), 

(V2:,V^), 

where c G { — 1, 0, 1} and a, 6 G M. Moreover, in the 
case c = we can scale the coefficient b to obtain 
b G {-1,0,1}. Additionally, if c = 6 = in the 
ninth subalgebra then we can set a G { — 1,0,1}. 

4 Group- invariant solutions 

Selected group-invariant solutions of (HI) have been 
studied in [TD] but without reference to classes of in- 
equivalent subalgebras and without noting that it 
is possible to set (3 = 0. Consequently, some of the 
ansatze presented in [TU] , which lead to a reduction 
of the number of independent variables of ([1]), are 
overly intricate. More precisely, they could be re- 
alised by means of considering reduction using one 
of the inequivalent subalgebras ([6|) or ^ and sub- 
sequently acting on the resulting invariant solutions 
by finite symmetry transformations. 

We investigate potentially interesting group- 
invariant reductions upon using ansatze that are 
based on the above classification of subalgebras. 
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If possible, we relate them to the solutions given 
in [TU]. 

4.1 Reductions with one-dimensional 
subalgebras 

Here wc give the complete list of reduced equations 
of H]) with the parameters F 0, (3 = based on 
the classification of inequivalent subalgebras ^ . In 
what follows V and w are functions of p and q. 

1. (D + aVr). Suitable invariants of this subalgebra 
for reduction are 

p— a;cos(alni) + ?/sin(aln<), 
q = — a;sin(alnt) + ?/cos(alni), 
V = tip. 

Using them as new variables, the vorticity equation 
dl]) is reduced to: 

w - a{qwp - pwq) - F{v - a{qVp - pVq)) + 

VqWp - VpWq = 0, 
W = Vpp + Vqq. 

2. {T> + avj.). Invariants of this subalgebra are 

p = X — ahit^ 

q = y, 

V = tip, 
which reduce ([T|) to 

W + aWp — F{v + aVp) ~ VpWq + VqWp = 0, 
W = Vpp + Vqq. 

3. (v-r ± Vf + av.^). Defining e — ±1, wc find the 
following invariants: 

p~ X coset + y sinet, 
q = —X sin et + y cos st, 
V ^ — eat. 

In these variables, ^ reads 

£{qwp - pWq) - sF{qVp - pVq + a) + 

VpWq - VqWp = 0, 
W = Vpp + Vqq. 

4. (vr + cv^). Here we have the invariants 

p = v^.t2 + 
q = t, 

X 

V = tp + c arctan — . 



This gives the same ansatz as was used in [T^ (Case 
2) upon additionally applying transformation ((31). 
The corresponding reduced equation reads 

c 

w„ — Fvn w„ = 

p 

1 

W = Vpp + -Vp. 

p 

5. (vt + + cv^). Invariants of this subalgebra 
are 

p ^ X — at, 

V = ip — ct. 

This is an ansatz for a traveling wave solution in 
x-direction, which includes the well-known Rossby 
waves. In |10| a similar ansatz was combined with 
a traveling wave ansatz also in y-direction (Case 3). 
However, as was indicated above, the additional 
consideration of waves in j/-direction is not neces- 
sary at this stage of reduction. The equation corre- 
sponding to the above ansatz is 

aWp — F{aVp — c) — VpWq + VqWp = 0, 
W = Vpp + Vqq. 

6. (vj; + cv^). A suitable ansatz for the invariants 
of this subalgebra is provided by 

p = x, 

q = y, 

V = ijj ~ ct. 
In this case ([!]) is reduced to 

Fc — VpWq + VqWp = 0, 
W = Vpp + Vqq. 

7. (v^). It is not possible to make an ansatz for tp 
for this subalgebra in the framework of the classical 
Lie approach. Hence, no reduction can be achieved 
through the gauging operator v^. 

4.2 Reductions with two-dimensional 
subalgebras 

To give also an example for a reduction using a two- 
dimensional subalgebra, let us consider the algebra 
(vj + avj. + cv^,, Vj^ + b'Vjp). In what follows v is a 
function of p. The invariants of this algebra are 

p = X — at, 

, , c (8) 

V = tp — by X, 

a 

provided that a 7^ 0. The corresponding reduced 
ODE of (P) then reads 

(a -I- b)vppp — Favp = 



4 



with the general solution 



V = Vi cxp 



Fa 



W2exp 



Fa 



a + b 



P] + V3, 



where Vi = const, i = 1,2,3. Transforming back 
to the original variables, renaming the constants Vi 
and applying the transformation ([¥]), we obtain the 
invariant solution 



V'l exp I \/ (x + ^t- at) j + 



which for Fa/ {a + b) < gives rise to a travelling 
wave solution. 

For the singular case a = 0, we cannot use ansatz 
([8|). Instead, we have the ansatz 

V = — ct ~ by 
and the reduced vorticity equation reads 



Fc + bVn 







which gives rise to a polynomial solution in the case 
7^ 0. If 6 = 0, we get the condition that c = and 
the ansatz ip = v{x) itself is the solution of ([T]). 

5 Summary and further 
comments 

In the sections above, we discussed distinct cases 
of reduction by using inequivalent Lie subalgebras. 
The main advantages of this systematic approach 
are the following: 

• Simplification of equation ([T]) since we only 
have to consider the case F ^ 0, P = Q. 

• Considering a minimal number of essential Lie 
subalgebras for reduction. 

• Simplifying the ansatz for the reduced equa- 
tions. 

• Optimal preparation of the reduced equations 
for further investigations. 

Upon using this approach, we have discussed all 
possible reductions by means of one-dimensional Lie 



subalgebras. We have to note that it is not possi- 
ble to use the subalgebra {v^) for obtaining group- 
invariant solutions, since in this case there is no way 
to make an ansatz for ijj. In principle, the remain- 
ing two-dimensional Lie subalgebras can be used for 
reduction as well. 

Moreover, the differential equations obtained by 
reduction could again be investigated by means of 
symmetry techniques. In general, some of the sym- 
metries of these reduced equations will be induced 
by the symmetries of the original equation. How- 
ever, sometimes there may be additional symme- 
tries that are not induced in this way [12l [H] and 
which are called hidden symmetries [1]. They usu- 
ally play an important role in the study of differ- 
ential equations, as they may allow to reduce equa- 
tions further than initially expected. 

We may also note that it is still possible to gen- 
eralise some of the results of this paper. By con- 
sidering eqn. Il]) as a system of two PDEs in the 
two dependent variables tp and it is possible to 
construct partially invariant solutions [15j . For this 
class of exact solutions, one needs at least two de- 
pendent variables. For the first set of these depen- 
dent variables, it is possible to introduce new in- 
variant variables, for the second set we keep the 
old noninvariant variables. The resulting reduced 
equations then also split in two sets of equations 
which have to be solved one after another. For this 
purpose, we could e.g. use subalgebras containing 
the operator v^. In this case, it is still not possi- 
ble to make an ansatz for ip but it is possible to do 
so for (. The resulting reduced system of differen- 
tial equations may give rise to a much wider class 
of exact solutions than pure group-invariant solu- 
tions. An investigation of this class of solutions for 
the case F = 0, f3 ^ will be given elsewhere. 
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